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PREFACE 
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ratory  (CERL) .  Dr.  Ashok  Kumar  of  the  CERL  was  the  primary  point  of  contact 
for  the  work.  The  report  was  written  and  prepared  by  Dr.  Thomas  M.  Walski, 

Mr.  Wayne  W.  Sharp,  and  Dr.  F.  Douglas  Shields.  Jr.,  of  the  Water  Resources 
Engineering  Group  (WREG) ,  Environmental  Engineering  Division  (EED) ,  EL. 

Ms.  Cheryl  Lloyd  of  the  WREG  proofread  the  report. 

The  work  was  accomplished  under  the  direct  supervision  of  Dr.  Shields 
and  Dr.  Paul  R.  Schroeder.  who  served  as  Acting  Chief.  WREG,  and  Dr.  R.  L. 
Montgomery,  Chief,  EED,  and  under  the  general  supervision  of  Dr.  John  Keeley, 
Assistant  Chief,  EL,  and  Dr.  John  Harrison,  Chief,  EL. 

COL  Dwayne  G.  Lee,  CE,  was  the  Commander  and  Director  of  WES. 

Dr.  Robert  W.  Whalin  was  Technical  Director. 

This  report  should  be  cited  as  follows: 

Walski,  Thomas  M. ,  Sharp,  Wayne  W. ,  and  Shields,  F.  Douglas,  Jr. 

1988.  "Predicting  Internal  Roughness  in  Water  Mains,"  Miscellaneous 
Paper  EL-88-2,  US  Army  Engineer  Waterways  Experiment  Station, 

Vicksburg,  Miss. 
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PREDICTING  INTERNAL  ROUGHNESS  IN  WATER  MAINS 


PART  I:  INTRODUCTION 


Background 


1 .  Some  measurement  of  the  internal  roughness  of  pipes  is  an  important 
parameter  entering  into  virtually  all  calculations  Involving  the  sizing  and 
analysis  of  water  distribution  systems.  In  water  utility  practice  in  the 
United  States,  the  Hazen-Williams  C-factor  is  the  most  commonly  used  parameter 
to  represent  the  carrying  capacity  (and  internal  roughness)  of  water  mains. 

2.  In  modern  cement  mortar  lined  pipes  and  plastic  pipes,  the  pipe 
roughness  changes  very  slowly  over  the  life  of  a  pipe  (with  the  exception  of 
waters  with  significant  scaling  potential  or  poor  removal  of  aluminum  hydrox¬ 
ide  floes).  However,  there  exist  many  miles  of  old,  bare  metal  water  mains  in 
which  the  pipe  roughness  has  changed  and  continues  to  change  significantly 
over  time. 


Estimation  of  C-Factors 

3.  Knowledge  of  the  pipe  roughness  (or  C-factor)  of  in-plare  water 
mains  is  critical  for  pipe  sizing  calculations,  and  there  are  several  ways  of 
obtaining  this  information.  The  first  is  to  use  "typical  literature"  values 
for  pipe  roughness.  As  will  be  shown  later  in  this  paper,  the  only  thing 
"typical"  about  such  values  is  that  they  vary  widely  from  water  system  to 
water  system. 

A.  A  second  method  is  to  determine  the  internal  roughness  of  existing 
mains  by  calculating  the  value  of  pipe  roughness  that  will  make  a  computer 
model  of  the  pipe  network  appear  to  be  calibrated  over  a  wide  range  of  condi¬ 
tions.  Some  examples  of  such  a  procedure  are  given  by  Walski  (1986)  and 
Ormsbee  and  Wood  (1986).  These  methods  essentially  give  an  "effective" 
C-factor  because  many  small  pipes  and  minor  losses  are  eliminated  when  the 
distribution  system  is  simplified  (skeletonized)  to  make  it  more  workable. 


5*  A  third  method  Is  to  actually  measure  the  pipe  roughness  in  the 
field  by  conducting  head  loss  tests  (Walski  1984).  This  method  is  quite 
expensive  but  gives  the  most  credible  results. 

6.  None  of  the  above  methods  for  determining  C-f actors ,  however,  pro¬ 
vide  a  means  for  extrapolating  pipe  roughness  into  the  future.  This  paper 
will  review  pipe  roughness  data  collected  from  a  wide  variety  of  sources  over 
many  years  and  provide  a  theoretically  based  procedure  for  predicting  future 
C-factor  values.  The  procedure  will  require  information  on  historical 
C-factors  or  water  quality. 


PART  II:  SOLUTION  OF  HEAD  LOSS  EQUATIONS  FOR  C 


Use  of  Hazen-Wllliams  Equation  for  Rough  Flow 


7.  The  Hazen-Williams  equation,  which  Is  traditionally  used  In  the 
United  States  for  calculating  head  loss  and  flow,  Is,  strictly  speaking,  only 
applicable  to  smooth  flow  (i.e.,  flow  in  which  roughness  elements  in  the  pipe 
wall  do  not  penetrate  through  the  laminar  bourdary  layer) .  The  equation 
becomes  somewhat  Inaccurate  for  transition  flow  and  even  less  accurate  for 
completely  rough  flow  which  is  usually  the  case  in  old  unlined  water  mains. 
The  equation  can  be  written  as: 


V 


0.55  C  D 


0.63 


s0.54 


(1) 


where 

V  ■  velocity,  ft/sec 
C  •  Hazen-Williams,  C-factor 
D  -  diameter,  ft 
S  »  hydraulic  gradient,  ft/ft 

8.  The  Hazen-Wllllams  equation,  however,  is  still  used  even  in  rough 
flow  because  the  error  in  predicting  head  loss  is  not  significant  except  for 
long  pipes  with  very  high  velocity.  In  most  water  distribution  system  prob¬ 
lems,  the  Hazen-Wllliams  equation  is  sufficiently  accurate. 

9.  In  rough  flow  the  C-factor  becomes  a  function  of  the  velocity  and 
can  be  corrected  by  the  equation  (Walski  1984)  given  below 

C  -  C  (V  /V)0,081  (2) 

o  o 


where 

C  •  C-factor  at  velocity  V 

C  «  measured  C-factor  at  velocity  V 
o  o 

V  *  actual  velocity,  L/T 
Vq  -  velocity  at  which  C-factor  measured,  L/T 
The  velocity  at  which  the  C-factor  is  measured  can  vary  widely,  but  a  typical 
value  of  roughly  3  ft/sec  (0.9  m/sec)  can  be  used  for  most  calculations. 
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Equation  2  shows  that  the  error  associated  with  assuming  a  constant  but 
typical  velocity  is  minor. 


Use  of  Darcy-Weisbach  Equation  for  Head  Loss 

19.  The  more  theoretically  correct  equation  for  head  loss  is  the  Darcy- 
Welsbach  equation  which  can  be  solved  for  velocity  as 


V  - 


(3) 


where 

2 

g  -  acceleration  due  to  gravity,  ft/sec 
f  -  Darcy-Weisbach  friction  factor,  dimensionless 
The  friction  factor  can  be  calculated  based  on  properties  of  the  flow  and  the 
wall  roughness. 


Relationship  Between  Roughness  and  C-Factor 

11.  Because  of  the  fact  that  the  C-f actor  is  a  function  of  velocity, 
the  initial  discueslon  below  will  center  on  pipe  roughness  (e) ,  which  repre¬ 
sents  the  height  of  equivalent  sand  grain  roughness.  When  C-f actors  are  dis¬ 
cussed,  they  will  be  in  regard  to  some  known  velocity  (say  3  ft/sec 

(0.9  m/sec)) . 

12.  By  solving  both  the  Hazen-Wllllams  and  Darcy-Weisbach  Equations  for 
head  loss  and  eliminating  terms,  the  C-f actor  can  be  related  to  the  Darcy- 
Weisbach  friction  factor  by 


C 


17.25 


f°-5*(VD)0-081 


(A) 


The  friction  factor,  f  ,  can  be  related  to  internal  pipe  roughness  over  a 
wide  range  of  flows  by  the  Colebrook-White  transition  equation  (Colebrook 
1939).  However,  this  equation  cannot  be  solved  explicitly  for  f  .  Since 
this  paper  concerns  pipes  with  significant  tuberculatlon  and  scale,  a  simpli¬ 
fied  form  of  the  Colebrook-White  equation,  attributed  to  von  Karman  and 


mfwuwwviniuwiwiutNi 


Prandtl  (von  Karman  1930) ,  that  la  applicabla  only  to  fully  rough  flow,  can  ba 
used.  This  equation  can  be  solved  explicitly  for  f  as  given  below. 

f  -  [1.14  -  2  log  (e/D)]“2  (5) 

where  e  is  the  roughness  height  (same  units  as  D) .  (All  logarithms  In  this 
paper  are  base  10.)  Equation  5  can  be  substituted  into  Equation  4  to  give  a 
formula  relating  roughness  height  to  C-factor  for  a  velocity  of  3  ft/sec 
(0.9  m/sec)  as  shown  below  (neglecting  the  D8*®8*  tenn) 

C  -  [14.6  -  25.6  log  (e/D)]1,08  (6) 

Equation  6  is  used  below  to  make  conversions  between  C-factors  and  roughness 
height . 

13.  Equation  6  appears  as  a  nearly  straight  line  on  a  semi-log  plot 
over  a  typical  range  of  values  for  e/D  .  The  exponent,  1.08,  makes  the  equa¬ 
tion  somewhat  difficult  to  use  and  can  be  eliminated  by  determining  the 
straight  line  that  most  closely  approximates  Equation  6  over  the  range  of 
C-factors  of  Interest.  This  line  can  be  given  by 

C  -  18.0  -  37.2  log  (e/D)  (7) 

Equation  7  will  be  used  in  the  methods  presented  later  to  predict  C-factors  as 
a  function  of  age  of  pipes  and  water  quality. 

14.  While  the  complete  Colebrook-White  transition  formula  cannot  be 
solved  explicitly  for  f  over  a  wide  range  of  Reynolds  numbers  and  pipe 
roughness,  other  formulas,  most  notably  one  by  Swamee  and  Jain  (1976)  can  be 
solved  explicitly  for  friction  factor  for  smooth,  transition  and  rough  flow. 
The  formula  of  Swamee  and  Jain  is  given  below 


f  - 


0.25 


0. 27e  . 

~d“  + 


(8) 


where  N  is  the  Reynolds  number. 
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15.  Equation  8  can  ba  lnaartad  into  Equation  4,  and  using  a  velocity  of 
3  ft/sac  (0.9  m/sec)  and  eliminating  the  Reynolds  number  term  which  is  negli¬ 
gible  for  rough  flow,  the  C-factor  can  be  given  by 

C  -  33.3  [log  (e/3.7D)2]0,54  (9) 

The  above  equation  can  be  simplified  by  realising  that  the  log  term  (even 
though  it  is  negative)  is  always  raised  to  the  second  power  so  that  the  abso¬ 
lute  value  of  the  log  can  be  taken  and  the  exponents  can  be  combined  to  give 

C  -  33.3  |  log  (e/3. 7D)  |  1,08  (10) 

Equations  6,  7,  9,  and  10  give  virtually  the  same  value  for  C-factor  for  any 
value  of  pipe  roughness  greater  than  0.001  ft  (0.305  mm)  and  are  reasonably 
close  for  values  to  0.0001  ft  (0.0305  nan). 


PART  III:  PREDICTING  GROWTH  OP  ROUGHNESS  IN  PIPES 


Early  Observations 

16.  Tha  growth  of  pipe  roughness  with  tine  has  been  observed  for  some 
time.  Williams  and  Hasan  (1920)  in  their  original  hydraulic  tables  give  sug¬ 
gested  values  for  C-factor  for  each  diameter  pipe  for  various  ages.  They  do 
not  quantitatively  describe  the  effect  of  water  quality  on  pipe  roughness. 
Nevertheless,  their  suggested  values  have  been  used  as  virtually  "gospel"  by 
engineers  since  they  were  first  published. 

17.  The  work  of  von  Karmen  (1930),  Prandtl  (1933),  and  Nlkuradse  (1932) 
showed  that  friction  factors  in  the  Darcy-Weiebach  equation  could  be  related 
to  the  pipe  wall  roughness  and  the  Reynolds  number  of  the  flow.  Colebrook 
(1938)  developed  an  analytical  relationship  to  describe  this. 

18.  More  importantly  for  this  paper,  Colebrook  and  White  (1937) 
addressed  the  problem  of  changing  pipe  roughness  with  time.  Using  their  data 
and  data  prepared  by  a  Committee  of  the  New  England  Water  Works  Association 
(1935),  Colebrook  and  White  advanced  the  hypothesis  that  pipe  roughness  grows 
roughly  linearly  with  time  and  that  the  rate  of  growth  depends  most  highly  on 
the  pH  of  the  water.  They  assumed  that  roughness  was  virtually  aero  when  the 
pipe  was  new.  A  slightly  modified  form  of  their  equation  to  predict  roughness 
is  given  in  Streeter  (1971)  as 


e  ■  e  +  at 

o 


(ID 


where 

e  -  absolute  roughness  height,  L 
eQ  "  roughness  height  at  time  aero,  L 
a  ■  growth  rate  in  roughness  height,  L/T 
t  ■  time,  T 

19.  The  hypothesis  that  e  is  a  linear  function  of  t  is  examined  in 
detail  below.  Equation  11  can  now  be  substituted  into  Equations  6,  7,  9,  and 
10  to  give  a  method  of  predicting  C-factors  as  a  function  of  time.  This  is 
done  for  Equations  7  and  10  below: 
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Colebrook-Whlte 


C  -  18.0  -  37.2  log  (X) 


(12) 


Swamee-Jain 


C  -  33.3  |  log  (0.27X)  |  1,08 


(13) 


where  X  »  (e  +  at.)/D  . 

0 

20.  Equations  12  and  13  contain  two  constants.  eQ  ,  tha  Initial  rough¬ 
ness.  and  a  ,  tha  roughness  growth  rata,  which  must  ba  determined  empiri¬ 
cally.  The  initial  roughness  depends  on  the  pipe  material,  but  a  typical 
value  of  0.0006  ft  (0.18  mm)  gives  reasonable  results  for  new  metal  plpaa  and 
is  reasonably  close  to  values  reported  by  Lamont  (1981).  Methods  to  predict 
the  roughness  growth  rate  must  take  into  account  water  quality.  Data  on  the 
relationship  of  roughness  growth  to  water  quality  are  presented  below. 

Historical  Data  on  Roughness  Growth  Rate 

21.  Colebrook  and  White  (1937)  reported  values  of  roughness  growth  rate 
ranging  from  0.000018  ft/yr  (0,066  iran/yr)  to  0.00017  ft/yr  (0.63  mm/yr). 

Based  on  data  they  obtained  from  the  New  England  Water  Works  Association 
(1935),  they  proposed  the  following  equation  to  predict  the  growth  rate, 
although  they  admitted  chat  it  was  "little  better  than  a  guess." 

a  -  0.0833  exp  (1.9  -  0.5  pH)  (14) 

22.  Another  significant  finding  of  Colebrook  and  White  was  that  the 
loss  of  carrying  capacity  was  due  much  more  to  the  increase  in  pipe  roughness 
rather  than  the  decrease  in  pipe  diameter  due  to  the  space  occupied  by  the 
roughness  elements. 

23.  The  roughness  growth  rates  that  could  be  back  calculated  from  the 
tables  of  Williams  and  Hazen  (1920)  varied  considerably  depending  on  the  diam¬ 
eter  but  generally  corresponded  to  0.002  ft/yr  (0.6  mm/yr).  Their  values  were 
based  on  a  fairly  limited  number  of  observations  and  they  noted,  "In  general 
it  may  be  stated  that  rather  large  deviations  from  the  indicated  rates  of 
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reduction  in  carrying  capacity  are  found  In  Individual  cases,  but  that.  In  the 
experience  of  the  authors,  the  variations  are  about  as  often  in  one  direction 
as  in  the  other o"  They  recommend  field  head  loss  tests  when  a  high  degree  of 
accuracy  is  required. 


24.  As  the  years  went  by,  more  investigators  reported  data  on  growth  of 
internal  pipe  roughness  versus  time.  Lamont  (1981)  presented  the  most  thor¬ 
ough  compendium  of  data  on  pipe  roughness  and  gave  the  most  rational  guidance 
on  the  effects  of  water  quality  on  roughness  growth.  His  approach  consisted 
of  identifying  four  "trends"  in  roughness  growth  and  determined  a  linear 
growth  rate  for  each  trend  as  given  in  Table  1 .  Lamont  used  the  Langelier 
Index,  which  is  essentially  an  indicator  of  the  saturation  of  the  water  with 
respect  to  calcium  carbonate,  to  characterize  water  quality,  numerical  values 
for  growth  rate  as  compared  with  Langelier  Index  are  shown  in  Table  1  and  can 
be  predicted  by  Equation  13  below 


a  -  10 


-(4.08+0.38  LI) 


for  LI  <  0 


where  LI  is  the  Langelier  Index.  Note  that  Equation  15  is  developed  only 
for  Langelier  Indexes  less  than  zero  (i.e.  corrosive  water)  and  svould  not  be 


Table  1 


Roughness  Growth  Rate  for  Varying  Water  Qualit 


Growth  rate 

tnm/yr 

Langelier 

Trend 

Name 

(ft/yr) 

Index 

1 

Slight  attack 

0.025 

0.000082 

0.0 

2 

Moderate  attack 

0.076 

0.00025 

-1.3 

3 

Appreciai.  '  e  attack 

0.25 

0.00082 

-2.6 

4 


Severe  attack 


0.76 

0.0025 


-3.9 


extrapolated  to  waters  with  scale  forming  tendencies  (l.e.  positive  Langelier 
Index)  or  waters  with  other  postprecipitation  problems.  The  range  of  values 
for  growth  rate  of  roughness  reported  by  Lamont  Is  similar  to  that  reported  by 
Colebrook  and  White. 

25.  While  it  is  possible  to  predict  the  roughness  growth  rate  for 
waters  with  a  negative  Langelier  Index,  there  is  considerably  less  data  avail¬ 
able  to  predict  this  rate  for  positive  Langelier  Index  values.  Natural  waters 
usually  do  not  have  a  positive  Langelier  Index,  and  therefore,  the  problems 
with  supersaturated  water  usually  only  exist  in  hot  water  lines  or  for  utili¬ 
ties  which  soften  water  but  do  not  recarbonate.  Larson  and  Sollo  (1967) 
describe  a  city  in  which  2  in.  (50.8  mm)  deposits  of  magnesium  hydroxide  were 
found  in  8  in.  (203.2  mm)  pipes  and  another  in  which  the  C-f actor  dropped  from 
135  to  100  in  15  years  due  to  calcium  carbonate  and  magnesium  hydroxide  depo¬ 
sits.  They  were  able  to  reproduce  this  phenomenon  in  stainless  steel  pipes  in 
the  laboratory. 

26.  Another  problem  is  post  precipitation  of  aluminum  hydroxide  in  dis¬ 
tribution  systems  where  the  water  is  treated  using  alum.  No  quantitative 
model  is  available  to  relate  loss  of  carrying  capacity  to  pH  or  aluminum  con¬ 
centration,  but  loss  of  carrying  capacity  has  been  observed  by  Costello 
(1982),  Emery  (1980)  and  Larson  and  Sollo  (1967).  Emery  reported  that 
C-factors  dropped  from  140  to  a  range  of  109  to  100  over  a  10-  to  23-year 
period  due  to  such  deposits.  This  corresponds  to  a  roughness  growth  rate  of 
0.00025  ft/yr  (0.076  mm/yr) . 

27.  A  committee  on  Loss  of  Carrying  Capacity  of  Water  Mains  of  the 
California  Section  American  Water  Works  Association  (1962)  presented  data  on 
changes  in  C-factor  for  the  Los  Angeles  water  system.  They  presented  the 
results  of  over  70  tests  for  pipes  with  diameters  between  4  and  16  in.  (101.6 
and  406.4  mm).  Roughness  growth  rates  were  calculated  for  each  diameter  and 
were  found  to  range  from  0.0020  to  0.0011  ft/yr  (0.61  to  0.34  mm/yr). 

28.  Hudson  (1966,  1966,  1973)  reported  the  results  of  numerous  head 
loss  tests  in  nine  cities.  While  the  data  were  reported  in  terms  of  C-factor 
instead  of  roughness,  and  there  was  a  considerable  amount  of  scatter  among  the 
data,  it  was  still  possible  to  back  calculate  roughness  growth  rates,  using 
Equation  6  for  each  city  as  shown  in  Table  2. 

29.  The  above  data  show  that  most  surface  waters  tend  to  be  corrosive. 
However,  the  limestone  well  water  in  San  Antonio  proved  to  be  much  less 
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Table  2 

Roughness  Growth  Rates  Calculated  from  Hudson's  Data 


City 

Growth  Rate 
mm/yr 
(ft/yr) 

Description  of  Water  Quality 

Atlanta 

0.61 

Soft  river  water 

0.0020 

Fort  Worth 

0.55 

Carrying  capacity  lost  even  in  concrete 

0.0018 

and  cement  mortar  lined  mains 

Denver 

0.18 

Mountain  reservoirs 

0.0006 

New  Orleans 

0.16 

River  water 

0.00052 

Cincinnati 

0.14 

River  water 

0.00043 

Chicago 

0.10 

Lake  water,  alum  treated 

(south) 

0.00033 

St.  Paul 

0.045 

Unsoftened  surface  water 

0.00015 

Chicago 

0.027 

Lake  water,  no  alum 

(north) 

0.0009 

San  Antonio 

0.015 

Wells  in  limestone 

0.00005 

corrosive  because  it  Is  near  saturation  with  respect  to  calcium  carbonate. 
Unfortunately,  data  on  the  pH,  calcium,  and  alkalinity  of  the  waters  in  the 
above  table  were  not  available  so  a  quantitative  method  to  predict  growth  rate 
as  a  function  of  water  quality  could  not  be  developed.  The  higher  growth  rate 
for  the  Chicago  south  treatment  plant  is  attributed  to  postprecipitation  of 
aluminum  hydroxide. 


Examination  of  Linear  Hypothesis 

30.  If  sequential  measurements  of  C-f actor  spanning  several  years  are 
available  for  the  pipe  sizes  and  water  source  of  interest,  and  if  water  qual¬ 
ity  has  been  fairly  unchanged  over  this  period  of  record,  then  the  engineer 
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can  predict  C-factors  as  a  function  of  pipe  age  using  Equation  12  or  13.  The 
likely  error  Inherent  in  such  an  approach  is  probably  acceptable  if  the  hypo¬ 
thesis  that  pipe  roughness  is  a  linear  function  of  pipe  age. 

31.  The  predictive  equations  for  C  (Equations  12  and  13)  are  based  on 
the  hypothesis  that  pipe  roughness,  e  ,  is  a  linear  function  of  pipe  age, 

t  .  Tills  hypothesis  was  examined  using  all  available  data.  The  data  set  used 
consisted  of  sequential  measurements  of  C  for  pipes  of  known  diameter  and 
age.  Some  319  data  points  from  seven  utilities  plus  values  from  Williams  and 
Hazen  (1920)  and  average  values  from  Lamont  (1981)  were  used.  First,  all 
C-factors  were  converted  to  pipe  roughness  values  (e)  using  Equation  16. 

Then,  for  each  utility,  pipe  roughness  growth  rates  (a)  were  determined  for 
each  pipe  diameter  by  linear  regression.  The  values  of  e  computed  with 
Equation  16  were  input  to  the  regression  as  dependent  variable  values,  and 
pipe  ages  (t)  were  the  dependent  variable  values.  The  a  values  determined 
from  regression  were  then  combined  with  published  or  measured  values  of  eQ 
and  the  known  values  of  D  and  t  to  compute  X  . 

32.  The  amount  of  scatter  present  in  the  entire  data  set,  and  thus  the 
validity  of  the  linear  hypothesis,  was  examined  by  plotting  C  versus 

log  (X)  as  shown  in  Figure  1.  Figure  1  indicates  that  a  linear  function  of 
t  is  a  good  predictor  of  pipe  roughness  and  therefore  Equations  12  or  13  may 
be  used  to  predict  C-factors  as  a  function  of  pipe  age.  The  95-percent  confi¬ 
dence  interval  about  the  regression  line  is  also  shown  in  Figure  1.  Predicted 
C-factors  are  within  an  interval  of  ±15  of  observed  values  with  95-percent 
confidence.  The  coefficient  of  determination  for  both  equations  is  0.87, 
which  means  that  the  equations  describe  87  percent  of  the  variation  in  C. 

This  is  especially  good  considering  the  fact  that  the  data  were  collected  by  a 
wide  variety  of  investigators  for  a  wide  range  of  conditions.  For  example, 
different  investigators  use  different  methods  for  determining  minor  head 
losses. 

33.  In  addition  to  the  fact  that  all  of  the  observations  fall  reason¬ 
ably  close  to  both  predicted  curves,  the  two  equations  (12  and  13)  give  nearly 
equal  results  over  the  range  or  C  normally  encountered.  Therefore,  either 
can  be  used,  depending  on  the  preference  of  the  engineer.  Examples  of  their 
use  follow. 
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18  -  37.2  log(X) 


PART  IV:  APPLICATION  OF  PREDICTIVE  EQUATIONS 


Procedure 


34.  To  apply  the  predictive  equation  to  a  water  system,  the  engineer 
needs  to  first  determine  the  growth  rate  of  roughness.  This  can  be  done  in 
two  ways.  First,  if  C-factor  and  pipe  age  data  are  available,  the  engineer 
can  plot  pipe  age  versus  roughness  (e) .  Roughness  (e)  can  be  calculated  from 
C-factor  by  solving  Equation  6  or  10  for  e  : 


.-d(io[<i8-c)/37-21) 


(16) 


1.7  (d 


10-[CC°-926)/25.7j) 


(17) 


The  slope  of  the  e  versus  t  plot  will  be  the  growth  rate  which  can  be  used 
for  predicting  future  C-f actors.  The  intercept  with  the  vertical  axis  gives 
an  indication  of  initial  pipe  roughness.  However,  since  new  pipe  C-f actors 
usually  are  not  measured,  since  there  is  some  error  In  C-factor  measurements, 
and  since  the  predictive  equations  are  not  very  accurate  for  roughness 
approaching  zero,  the  estimates  of  Initial  roughness  may  appear  unreasonable 
in  some  cases  (e.g.  negative).  In  such  cases,  the  engineer  should  use  typical 
initial  roughness  values  of  approximately  0.0006  ft  (0.18  mm).  In  any  case, 
the  value  of  initial  roughness  will  usually  have  very  little  effect  on 
C-factors  after  the  first  few  years  of  a  pipe's  life. 

35.  If  the  transformed  data  do  not  fall  on  a  straight  line  with  reason¬ 
able  error,  the  engineer  must  consider  why  this  has  occurred.  The  usual 
explanation  is  that  water  quality  has  changed  during  the  life  of  some  of  the 
pipes.  The  engineer  must  then  try  to  determine  the  aging  rate  given  the  cur¬ 
rent  water  quality  since  that  water  quality  will  influence  future  changes  in 
C-factor. 

36.  If  historical  C-factor  values  are  not  available  or  water  quality  is 
expected  to  change  in  the  future,  then  the  engineer  must  predict  roughness 
growth  based  on  water  quality.  Equation  15  should  be  used  for  this  purpose 
since  the  Langelier  Index  is  based  on  more  than  simply  pH,  as  given  in  Equa¬ 
tion  14.  Corrosion  inhibitors  will  also  have  an  effect  on  growth  rate. 
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37.  Once  the  roughness  growth  rate  is  known,  the  C-factor  in  a  given 
year  can  be  predicted  using  Equations  12  and  13,  knowing  the  age  of  the  pipe 
for  that  year.  This  is  an  acceptable  procedure  if  the  precise  age  of  instal¬ 
lation  of  each  pipe  segment  is  known.  Determining  the  year  laid  for  each  pipe 
segment  can  be  a  fairly  tedious  process  and  is  one  that  is  currently  not 
required  in  most  pipe  network  models.  An  alternative  procedure  is  to  extrapo¬ 
late  the  C-factor  given  the  current  value.  The  current  value  can  be  deter¬ 
mined  from  field  tests  or  the  value  needed  to  calibrate  a  pipe  network  model 
as  discussed  earlier. 

38.  To  modify  the  predictive  equations  so  that  the  year  laid  need  not 
be  known,  it  is  necessary  to  first  substitute  for  pipe  age  in  the  equations  as 
shown  below 


t 


(18) 


where 


T  *=  year  of  interest  (e.g.  1995) 

Tq  -  year  of  installation 

It  is  possible  to  reorganize  Equations  12  and  13  to  eliminate  Tq  ,  which 
results  in  an  equation  to  predict  C-factor  in  any  year  T  ,  given  the  C-factor 
in  some  other  year,  T^  : 


C  =  18.0  -  37.2  log 


aCT-Tj,)  +  D  10(18“C1)/37"2 


(19) 


C  -  33.3 


I  log 


-C 


0.27a  (T  -  Tj)  +  D  10 


0.926/25.7 

1 


1.08 


(20) 


where 

C^  ■  C-factor  in  year  T^  (known) 

Tj  ”  year  in  which  C-factor  is  known 
T  -  year  in  which  C-factor  is  predicted 

39.  The  above  equations  can  be  easily  incorporated  into  computer  models 
of  pipe  networks  so  that  an  engineer  can  simulate  future  system  behavior. 


17 


Example  Problems 


40.  Two  example  problems  to  Illustrate  the  predictive  equations  are 
presented  below.  The  first  covers  the  case  In  which  no  C-f actor  data  are 
available,  but  the  year  the  pipe  was  laid  is  known.  The  second  covers  the 
case  in  which  there  are  some  historical  C-f actor  data. 

Example  1 

41.  It  is  necessary  to  predict  the  C-factors  in  1995  for  the  pipes 
described  in  the  first  three  columns  of  Table  3.  The  pipes  are  unlined  cast 
iron  and  Langelier  Index  for  the  water  is  -1.5.  Determine  C-factors  for  sev¬ 
eral  individual  pipes  plus  plots  of  C-factor  versus  age  for  6,  12,  24  and 

."c  in.  (152.4,  304.8,  635.0  and  914.4  mm)  pipes. 


Table  3 

Data  for  Example  1 


Pipe 

Diameter 

(ft) 

Year  Laid 

C  in  1995 

1 

0.5 

1925 

68 

2 

0.5 

1905 

64 

3 

0.5 

1940 

72 

4 

2.0 

1925 

91 

42.  The  rate  of  aging  can  be  estimated  from  Equation  15  as 

a  -  10-r4.08+0.38(-1.5)]  _  0>00031  ft/yr  (0.095  mm/yr) 

The  age  of  each  pipe  in  1995  can  be  determined  by  subtracting  the  year  laid 
from  1995.  Then,  using  an  initial  roughness  of  0.0006  ft  (0.18  mm).  Equa¬ 
tion  12  can  be  used  to  generate  the  C-factor  in  1995  as 

C  -  18.0  -  37.2  log  [0.0006  +  0.00031  (1995  -  T  )]/d| 
i  o  I 

A  plot  of  C-factors  versus  age  for  several  diameter  pipes  is  shown  in  Fig¬ 
ure  2.  Note  that  C-factor  decreases  more  rapidly  in  small  diameter  ripes. 
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Figure  2.  C-f actor  versus  age 


Example  2 

43.  In  this  problem,  several  C-factor  tests  were  conducted  and  the 
results  are  shown  in  the  first  three  columns  of  Table  4.  The  data  were  then 
transformed  as  shown  in  Equation  17,  and  the  res  : ts  are  shown  in  Figure  3. 

The  roughness  growth  rate  was  found  to  be  0.0019  ft/yr  (0.58  mm/yr)  using  lin¬ 
ear  regression.  Equation  20  was  then  used  to  predict  the  C-factors  in  2005: 


33.3 


1  log 

—  (C^*^^/25  7) 
0.0103  +  D  10  1  1  m*') 

1 

D 

1 

1.08 


Rate  of  Change  of  C-Factor 


44.  While  pipe  roughness  increases  linearly  with  time,  C-factors  ini¬ 
tially  decrease  rapidly  and  then  change  more  gradually  as  shown  in  Figure  2. 


0.0019 


1 _ I _ I _ L 

10  20  30  40  50 

t.  YR 

Figure  3.  Growth  rate  versus  time  for  Example  2 
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Table  4 

Data  for  Example  2 


Diameter 

C-factor 

Age  in 

e 

C  factor 

Plp« 

_  (ft) 

in  1985 

1985 

(ft) 

in  2005 

1 

2.0 

80 

20 

0.041 

69 

2 

1.0 

61 

35 

0.065 

54 

3 

1.0 

55 

45 

0.095 

50 

4 

0.5 

48 

45 

0.073 

41 

5 

0.5 

42 

50 

0.107 

37 

This  explains  why  results  of  pipe  cleaning  jobs,  when  not  accompanied  by 
cement  mortar  lining  or  a  change  in  water  quality,  are  usually  short  lived  as 
described  by  Dutting  (1968)  and  Frank  and  Perkins  (1955).  To  quantify  this 
effect,  it  is  possible  to  take  the  derivative  of  Equation  12  with  respect  to 
time.  This  gives  the  rate  of  change  in  C-factor  with  time. 


dC  -16.1 

dt  ’  [(eQ/a)  +  t)J 


(21) 


The  absolute  value  of  the  right  side  of  Equation  21  is  plotted  in  Figure  4  for 
roughness  growth  rates,  a  ,  of  0.001  and  0.0001  ft/yr  (0.305  and 
0.0305  mm/yr).  Figure  4  shows  th-t  the  C-factor  drops  off  initially, 
decreases  rapidly  and  then  changes  very  slowly. 
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PART  V:  SUMMARY  AND  CONCLUSIONS 


45.  The  aquae Ions  prasancad  in  this  paper  provide  an  easy-to-use  but 
accurate  method  for  predicting  C-factors  In  unllned  metal  pipes.  The  data 
substantiate  the  claim  of  Colebrook  and  White  (1937)  that  roughness  height  in 
such  pipes  grows  roughly  linearly  over  time. 

46.  The  equations  presented  depend  most  heavily  on  the  parameter,  a  , 
that  describes  the  rate  of  growth  of  roughness  over  time  and  depends  on  the 
quality  of  the  water  with  respect  to  corrosion  or  precipitation.  This  param¬ 
eter  should  be  determined  based  on  historical  data  for  the  given  water  system. 
If  this  information  is  not  available,  it  can  be  estimated  based  on  the 
Langeller  Index  or  extrapolated  from  other  utilities  with  similar  water 
quality. 

47.  Since  the  C-factor  is  a  logarithmic  function  of  time,  the  C-factor 
tends  to  decrease  most  rapidly  for  new  or  recently  cleaned  pipes.  However, 
after  a  few  decades  the  drop  in  C-factor  with  time  is  very  slow. 
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